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ABSTRACT 


>  A  dual-oanc  comparison  in  the  3WIH  and  LMIR  atmospheric  windows 
was  made  of  the  spatial,  temporal,  and  radiometric  statistics  of  clouds 
as  cbserved  Ground  Eased  Measurements  (GBM)  Sensor  of  the  Army 
Optical  Station  'AO^.  The  existence  of  the  GBM  archives  has  provided  a 
unique  opportunity  to  analyze  the  properties  of  cloud  images 
simultaneously  observed  in  the  3-5  and  the  3-12  micron  passoands.  It  is 
shown  that  the  level  of  correlation  in  clutter  between  the  two  passbands 


is  a  strong  function  of  the  meteorological  conditions.  Through 
observations  of  reentry  vehicles  passing  behind  clouds,  point 
measurements  of  the  attenuation  of  target  signatures  have  been  achieved. 
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1.  HRiooDcnoR 


1.1  TASK  SUMMARY 

This  report  is  submitted  to  Night  Vision  and  Electro-Optical 
Laboratories  as  the  final  report  for  Contract  Number  DAAB07-82-D-J009 , 
Subtask  157.  A  comparative  analysis  was  performed  to  determine  the 
level  of  correlation  in  radiometric  signature  that  exists  for  clouds 
simultaneously  observed  in  the  LWIR  and  the  SWIR  atmospheric  passbands. 
This  comparative  study  was  performed  on  nine  missions  of  the  GBM 
Archives  (Refer  to  Table  2-4  of  Section  2  of  this  report),  which  is  the 
cloud  data  processed  for  Task  I.  From  this  analysis,  it  has  been  shown 
that  in  these  cases  clouds  act  as  emitters  in  the  LWIR  and  as  scatterers 
in  the  SWIR.  It  was  also  shown  that  the  level  of  correlation  between 
passbands  is  a  strong  function  of  weather  conditions. 

A  summary  of  the  total  task  being  performed  by  Teledyne  Brown 
Engineering  (TBE)  in  this  effort  is  presented  here  to  provide  an  under¬ 
standing  of  the  intended  use  of  this  information.  This  effort  is  divid¬ 
ed  into  five  tasks,  and  these  tasks  have  subsequently  been  grouped  into 
two  phases.  Phase  I  (Tasks  I  and  IV)  was  to  determine  the  extent  of 
dual-band  cloud  data  accessible  from  the  GBM  archives  and  to  perform 
sufficient  analysis  on  these  data  to  show  the  level  of  correlation  (if 
any)  between  the  LWIR  and  SWIR  passbands.  Phase  II  (Task  II,  III,  and 
V)  will  be  a  continuation  of  this  analysis  on  additional  cases,  the 
development  of  a  technique  for  Che  conversion  of  these  data  into  a  form 
comparable  to  a  FLIR  image  and  ultimately  Che  development  of  a  wide- 
field-of-view,  dual-band  cloud  scene  based  on  this  statistical  analysis 
on  the  GBM  cloud  data. 

1.2  SYNOPSIS  OF  REPORT 

Section  2  of  this  report  gives  a  description  of  the  GBM  Data 
Base  and  the  extent  of  cloud  data  it  contains.  Since  several  sensors 
were  involved  in  data  collection,  they  are  also  described  in  this 


section. 


Section  3  is  a  detailed  discussion  of  the  theory/  of 
statistical  analysis  as  incorporated  in  this  study.  Included  in  this 
section  are  the  definitions  of  the  statistical  tools  used  in  this  study, 
and  a  discussion  of  the  choice  of  the  cosine  roll-off  apodization  i 

I 

function  applied  to  the  raw  data  stream  before  Fourier  analysis.  The 
final  choice  of  the  30-point  cosine  roll-off  is  shown  as  a  compromise 
between  high-frequency  aliasing  and  signal  amplitude  error.  Section  4 
contains  the  results  of  the  comparative  analysis  and  a  discussion  of  the 
potential  use  of  the  reentry  vehicle  as  a  calibrated  source  for  the 
determination  of  cloud  transmission  profiles.  Section  5  describes  the 
cloud  scene  generation  model  developed  by  TBE  for  Night  Vision  and 
Electro-Optics  Laboratories  and  suggests  a  method  for  the  temporal 
propagation  of  simulated  cloud  scenes. 

Appendix  A  of  this  report  describes  the  procedure  used  to  reduce 
Che  raw  PCM  data  and  to  convert  it  to  a  digital  form.  This  technique 
was  originally  created  to  reduce  reentry  vehicle  target  data  but  is  com¬ 
patible  with  the  reduction  of  cloud  data.  Appendix  B  discusses  Che  pro-  i 

cedure  for  calibration  of  the  GBM  sensor  during  mission  recording,  and 
describes  Che  measures  taken  Co  ensure  an  accurace  calibracion  standard. 

-Appendix  C  describes  the  extent  of  the  supporting  meteorological  data 
obtained  during  the  radiosonde  balloon  flights  at  the  time  of  each  33M 
mission  and  provides  temperature  and  dew  point  profiles  for  each  of  the  ^ 

missions  listed  in  Table  I.  This  appendix  also  describes  a  method  for-  I 

approximating  the  observed  cloud  ranges  based  on  these  profiles.  -Appen¬ 
dix  D  contains  polar  plots  of  the  boresighc  cracking  and  Sun  position  I 

for  the  above-mentioned  missions. 


2.  G8M  SENSOBS  AMD  DATA  BASE 


2.i  GBM  SENSOR  SYSTEM 

The  Ground  Based  Measurements  (GBM)  Sensor  System  is  a  passive, 
two-color  infrared  optical  scanning  and  tracking  sensor  mounted  on  a 
modified  NIKE-AJAX  azimuth-elevation  tracking  pedestal.  The  system  is 
located  at  the  Army  Optical  Station  (AOS)  on  Roi  Namur  (Figure  2-1)  in 
the  Kwajalein  Atoll.  The  system's  purpose  is  to  obtain  high-resolution 
infrared  signatures  and  other  temporal  and  spatial  characteristics  of 
reentry  targets  which  terminate  at  the  Kwajalein  impact  area. 

The  GBM  pointing  and  tracking  system  was  designed  to  operate  in 
any  of  several  modes: 

•  Radar  command 

•  TV  auto-track 

•  GBM  infrared  auto-track 

•  Manual  track 

•  Trajectory  retrace 

•  Console  control. 

These  modes  can  be  selected  by  means  of  a  track-mode  switch  on  the  Ope¬ 
rations  Console. 

The  sensors,  optics,  and  servo-drives  are  located  on  the  track¬ 
ing  pedestal  in  the  dome  of  the  AOS  building  shown.  The  Operations 
Console  and  other  electronics  are  located  in  a  34-ft  instrument  van. 

The  focal  plane  of  the  GBM  sensor  contains  two  staggered  line 
arrays  of  mercury-doped,  germanium  detectors,  which  are  cooled  bv  a 
closed-cycle,  gaseous-helium  system.  The  two  infrared  bandpasses  are 
obtained  by  filtering  each  array  to  a  different  bandpass.  The  system 
relative  response  for  each  waveband  is  provided  in  a  separate  document. 
Figure  2-2  shows  the  NEFD  of  each  GBM  passband  as  a  function  of  focal 
plane  irradiance.  There  are  31  detectors  for  LVIR  and  21  for  SWIR. 
Figure  2-3  displays  the  geometry  of  the  two  arrays.  The  sensor  employs 
a  dither  mirror  scanning  bidirectionally  at  a  rate  of  15  Hz.  Each 
detector  is  sampled  at  a  rate  of  25  kHz;  voltages  are  subsequently  digi¬ 
tized  and  encoded  for  recording.  Characteristics  of  the  GBM  sensor  are 


presented  in  Table  2-1.  The  GBM  sensor  focal  length  is  estimated  on  the 
basis  of  initial  measurements  made  ac  AOS. 


TABLE  2-1.  AOS  SENSOR  CHARACTERISTICS 


FRAME  RATE 
(frames/sec) 

FOCAL 

RATIO 

FOCAL  LENGTH 
(mm) 

— 

TIMING 

FIELD  OF  VIEW 
(mrad) 

GBM  Sensor 

30  scans 
FOV/sec 

f/6.9 

3J60 

IRIG  B 
IRIG  A 

2.2  by  5.3 

GBM  LLL  TV 

60 

f/2.8 

215 

GMT 

70,  diagonal 

35 -mm  Camera 

10 

f/3.7 

150 

IRIG  B 

120  by  168 

2.2  GBM  DATA  RECORDING 

The  GBM  system  records  data  using  analog  recorders,  video 
recorders,  and  film  cameras. 

2.2.1  Analog  Recorders 

The  digital  signal  outputs  for  each  color  band  are  independently 
recorded  on  two  Bell  and  Howell  lA-track  VR3700B  tape  recorders  with 
Bell  and  Howell's  enhanced  NRZ  electronics  for  the  PCM  channels.  The 
GBM  channel  assignment  is  shown  in  Table  2-II.  Channels  3  and  5  ara 
used  CO  record  the  prime  detector  data.  These  data  consist  of  sampled 
detector  outputs  that  have  been  converted  3-bic  digital  words.  Two  such 
data  streams,  and  hence  two  recorder  channels,  are  necessary  to  record 
all  of  the  detector  outputs  and  still  remain  below  the  bit  rate  capacity 
(4  Mbps  at  120  ips)  of  the  recorder. 

Both  PCM  data  channels  contain  a  series  of  29  eight-bit  .  data 
words,  totalling  '232  bits.  Three  of  the  29  words  are  sync  words.  The 
channel  3  sync  word  is  0000  0101  0000  11000  1101  1111.  The  channel  5 
sync  word  is  the  ones  complement  of  the  channel  3  sync  word.  The  first 
of  the  26  data  words  for  channel  3  is  the  output  from  detector  A4, 


DATA  DESCRIPTION 


.NME'i': 


RECORDING  MODE 


Edge  i 

Reserved  *or  Voice  A 

D'>ect 

T 

1 

Scanning  Secondary  Mirror  Position 

FM 

Analog  Voltage 

WB  Group  II 

2 

Spare 

TBD 

3 

GBM  Long-Wavelength  PCM 

PCM 

Radiometric 

Enhanced  NRZ 

4  - 

Spare 

TBD 

5 

GBM  Long-and  Short-Wavelength  PCM 

PCM 

Radiometric 

Enchanced  NRZ 

6 

Composite  Clock  (15  Hz  and  15  kHz) 

FM 

WB  Group  II 

7 

IRIG-B  Time  (Modified) 

FM 

WB  Group  II 

8 

Spare 

TBD 

9 

IRIG-A  Time  (Modi'*^ied) 

FM 

WB  Grouo  II 

10 

Spare 

TBD 

n 

GBM  Housekeeoing 

FM-PCM 

12 

Tape  Servo 

Direct, 

WB  Option  B 

13 

DC  Voltage  Levels,  Long-Wavelength 

FM 

WB  Group  II 

14 

DC  Voltage  Levels,  Short-Wavelength 

FM 

WB  Group  II 

Edge 

Voice  B 

Direct 

and  the  last  is  A29.  The  channel  5  data  words  actually  begin  with  the 
output  from  detectors  Al,  A2,  and  A3,  then  include  31  through  B21,  and 
end  with  A30  and  A31.  The  PCM  bit  rate  is  3.48  Mbps. 

A  time  reference  is  provided  by  FM-recording  IRIG-3  time  on 
channel  ~  and  IRIG-A  time  on  channel  9.  The  IRIG-3  time  code  contains 
100  pulses  per  second  with  a  carrier  frequency  of  1  kHz;  IRIG-A  contains 
1,000  pulses  per  second  with  a  carrier  frequency  of  10  kHz. 

The  housekeeping  data  are  FM-recorden  in  PCM  format  on  channel 
11  and  consist  of  44  eight-bit  words' at  a  bit  rate  of  approximately 
28.16  kbps.  Included  is  such  information  as  the  secondary  scanning  mir¬ 
ror  position,  pedestal  azimuth  and  elevation  position,  digitized  analog 
status  signals,  and  digital  status  signals.  In  addition,  there  is  a 
composite  15-Hz/15-kHz  clock  recorded  on  channel  6.  A  frequency  used  in 
tape  servo  compensation  is  direct-recorded  on  channel  12. 

2.2.2  Video  Recorders 

The  second  recording  subsystem  consists  of  three  IVC  Model  800 
video  recorders  which  record  the  display  from  the  GBM  LLL  TV  camera 
attached  to  the  GBM  sensor  telescope  mount  and  boresighted  with  the 
infrared  sensor,  a  scan-converted  display  of  both  infrared  bands  from 
the  GBM  sensor,  and  other  TV  signals  that  are  displayed  on  the  control 
console  monitors. 

The  video  recorders  serve  two  main  purposes.  The  first  is  to 
help  document  the  conditions  under  which  the  data  were  taken.  The 
second  is  to  aid  in  selecting  and  editing  the  data  to  be  formatted. 
Because  of  the  large  amount  of  recorded  material  obtained  during  a  mis¬ 
sion,  editing  is  necessary  to  minimize  the  amount  of  data  formatted  and 
to  time-correlate  significant  events. 

2.2.3  Film  Camera 

Besides  Che  video  recorders,  a  35-mm  motion  picture  camera  is 
boresighted  with  Che  sensor  for  additional  optical  recording.  Charac¬ 
teristics  of  Che  TV  and  35-mm  cameras  are  presented  in  Table  2-1. 


METSROLOGICAL  INSTRUMENTATION 


.  j 

The  GBM  Archives  includes  recordings  of  the  reentry  of  various 
strategic  warheads  and  their  associated  objects.  These  data  were  simul¬ 
taneously  recorded  with  several  measurement  instruments  including  radar. 
LLL  T/,  35-mm  camera,  and  the  GBM  dual-band  IR  sensor,  which  also 
provides  for  visual  interpretation  of  clouds.  In  addition  to  direct 
measurements  of  the  vehicle,  atmospheric  conditions  during  the  time  of 
reentry  were  measured  at  ground  level  and  with  radiosonde  balloon 
flights.  In  most  cases,  the  instrumentation  was  turned  on  for  test  runs 
(called  nominal  runs)  and  for  calibration.  During  these  periods,  cloud 
data  were  recorded  without  the  presence  of  the  reentry  vehicle  in  the 
FOV.  The  missions  cover  a  wide  variety  of  climatic  conditions  from 
clear  sky  to  heavy  overcast  and  raining.  The  missions  are  split 
approximately  equally  into  day  and  night  runs,  with  the  average  mission 
length  at  about  4  minutes. 

Table  2-III  gives  a  summary  of  the  extent  of  cloud  information 
and  the  respective  cloud  types  in  each  mission.  The  column  labeled  Scan 
Conv  refers  to  the  Scan  Converter  for  the  GBM  sensor  and  not  to  the  GBM 
sensor  data  itself.  Every  mission  listed  in  this  table  includes  the  GBM 
dual  IR  data  in  PCM  tape  format. 

Table  2-IV  shows  the  GBM  missions  as  they  appear  with  respect  to 
Sun  position:  the  boxed  case  numbers  are  those  chosen  for  further  analy¬ 
sis  and  subsequent  modeling  with  the  Extended  FOV  Cloud  Scene  Generator. 
These  missions  were  chosen  to  maximize  the  variation  in  time  of  day  and 
Sun  position,  and  they  cover  a  variety  of  climatic  conditions.  It  is 
noted  that  a  severe  limitation  of  this  data  set  for  ground-based  analy¬ 
sis  is  its  maritime  characteristics.  As  will  be  shown,  the  water  vapor 
content  significantly  reduced  the  SWIR  transmission  in  all  but  Case  60. 

The  designation  of  front-lit,  side-lit,  and  back-lit  clouds  is 
defined  as  shown  in  Figure  2-4,  in  which  Sun  rays  entering  within  60  deg 
of  the  sensor  boresight  produce  back-lit  clouds;  Sun  rays  entering 
between  60  and  120  deg  produce  side-lit  clouds,  and  120  to  180  deg 
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ARCHIVES 


CLOUD  INFORMATION:  (I)  Cl  LAN 

ALTITUDE:  (2)  HIGH  CIUUDS,  (J)  MIDltlM  CIOUDS.  (4)  LOW  CIOUDS 
DENSITY:  (5)  HEAVY  OVIRCASI,  (h)  MIDllIM  OVERCAST,  (7)  LIGHT  OVERCAST 
TYPE:  (8)  CIRRUS,  (9)  r.IJMIJIIIS,  (10)  SIKAIUS,  (11)  BROKEN  CLOUDS 
(12)  CLOUD  TYPE,  At  IIIOOI.  AND  (OVERAGE  WAS  NOT  DETERMINED 


CLOUD  INFORMATION:  (I)  CHAM 

ALTITUDE:  (2)  HIGH  CIOUDS,  (  J)  MEDIUM  Cl  DUDS,  (4)  LOW  CLOUDS 
DENSITY:  .  (5)  HEAVY  OVEMCASI,  (6)  MtDIUM  OVIRCAST,  (7)  LIGHT  OVERCAST 
lYPE:  (8)  CIRRUS,  (9)  CUMUIUS.  (10)  SIRAIUS,  (11)  BROKEN  CLOUDS 
(12)  CLOUD  TYPE,  AllllllDl,  AND  COVIKAGE  WAS  NOT  DETERMINED 


3.  THEORY  OP  STATISTICAL  ANALYSIS 


The  techniques  and  statistical  methods  used  to  perform  the  dual- 
band  analysis  in  this  report  are  defined  below.  Sections  3.1  through 
3.-*  describe  standard  statistical  methods  as  found  in  many  excellent 
references.  Section  3.5  provides  a  rationale  for  the  choice  of  the  30- 
point  cosine  roll-off  apodization  function,  used  on  all  data  before 
Fourier  analysis. 


3.1  IgAHS  AMD  VARIANCES 


While  not  of  primary  interest  in  a  comparative  analysis,  the 
determination  of  means  and  variances  of  the  input  data  is  required  for 
other  analysis  functions  such  as  cross-correlation.  For  a  series  of  N 
observations  (in  this  analysis  N=256),  the  most  probable  estimate  of  the 

mean  •_  is  given  bv , 

N 


y 


Xi 


i=l 


where  X,.  is  the  ith  data  element  of  the  X  sampled  data  points. 
The  oest  estimate  o*  the  variance  r”  is  liven  bv. 


N 


i=1 


where  N-1  is  usea  since  cnis  represents  the  number  of  degrees  of  freedom - 
left  after  determining  u  from  N  observations. 


3.2  CROSS-CORRELATIONS 

Given  two  mutually  stationary  N-channel  time  series,  the  cross¬ 
correlation  C  as  a  function  of  phase  lag  is  computed  by  the  formula, 

N  p 
t=j+1 


3-1 


where  k,  j=x«2, i  =  0,i,.--L;  L=  the  tnaximunj  number  of  lag 
pomes:  Xj  *-  is  channel  k  of  the  first  tine  series;  j'  13  channel  : 
of  the  second  time  series;  and  and  "y^  j  ^  are  the  corresponding 


3.3  SPECTRAL  DENSITIES  AND  CROSS-SPECTRA 

The  spectral  densities  S  and  cross-spectra  Cs  based  on  the  cor¬ 
responding  cross-correlation  function  are  calculated  from  the  formulas 

2^  r*(qAt)Cos  mq  At  aw  +  r*(mAt)Cos  aw 
q=l 

where  ra^®^  is  the  auto  correlation  function  at  zero  phase  lag,  At  is 
the  interval  between  samples,  and  r*  is  the  correlation  function  fil¬ 
tered  by  the  Bartlett  Method, 

r*  (kAt)  =  (i  -  r(kA t), 


CsU)  =^C(m)e-'^"’ 
m  =  (N-1) 

where  C'm)  represents  the  cross-correlation  between  data  sets  and  w 
the  spatial  or  temporal  frequenc^^  in  samples  per  data  interval. 


3.L  L£^/EL  OF  SIGNIFICANCE 

In  the  interpretation  of  cross-correlation  coefficients  and  sub¬ 
sequently  cross-spectra,  it  is  important  to  determine  whether  the  coef¬ 
ficient  values  are  likely  to  reflect  a  predictable  relationship  between 
data  sets.  A  useful  distribution  for  this  determination  is  the  proba¬ 
bility  PCr,N)  that  a  random  sample  of  N  uncorrelated  experimental  data 
points  would  yield  an  experimental  linear-correlation  coefficient  as 
large  as  or  larger  than  the  observed  value  of  jr|.  This  probability 
function  P(r,N)  can  be  computed  by 


“  ^-,-7i[(V-H)/2]  T  (-1)^  II 
^  r(v/2)  ■TMTrrr 


M  +  1 


for  even  values  of  v.  where  v  =  N  -  2  and  I  =  D.5  v  -  2'.  This  tech¬ 
nique  was  used  to  determine  the  levels  of  significance  discussed  in  Sec¬ 
tion  4. 


3.3  APODIZATION 

Because  of  Che  finite  nature  of  data  sets  and  the  manner  in 
which  Fourier  transform  routines  interpret  the  end  points  of  these  sets, 
it  is  necessary  to  smoothly  bring  the  end  points  to  zero  amplitude  if 
high-frequency  aliased  noise  is  to  be  avoided.  The  standard  technique 
used  to  smoothly  terminate  data  strings  is  to  multiply  them  by  an  apodi- 
z  at  ion  window.  The  form  of  the  window  chosen  for  this  study  was  the 
cosine  roll-  off  in  which  Che  M  samples  at  each  end  of  the  data  set  are 
multiplied  by  a  factor  from  one  to  zero  according  to  a  cosine  relation¬ 
ship.  The  following  discussion  gives  the  rationale  for  the  choice  of 
M  =  30. 


Figure  3-1  shows  the  Fourier  transform  of  a  6-cycle-per-data 
interval  sine  wave  having  been  filtered  by  a  cosine  roll-off  apodization 
window  of  the  indicated  width.  The  first  chart  of  this  figure  (M  =  0) 


corresponds  to  a  boxcar  window  or  no  apodization  beyond  starting  and 
stopping  the  data  string  to  define  the  sampled  data  set.  As  the  extent 
of  the  cosine  roll-off  increases,  the  high-frequency  lobes  or  feet  on 
the  transform  spike  begin  to  reduce,  but  as  this  occurs  the  aliased  dc 
noise  level  begins  to  increase  in  amplitude  and  spectral  width.  Fijure 
3-2  illustrates  the  increase  in  dc  noise  level  with  increasing  M.  Fig¬ 
ure  3-3  shows  the  same  relationship  of  the  percent  error  in  transform 
peak  amplitude.  Given  these  two  contributions  alone,  the  optimum  condi¬ 
tion  occurs  for  no  apodization  or  M  =  0.  However,  the  purpose  of  apodi- 
zation  is  to  remove  the  "feet"  produced  by  the  sudden  termination  of  a 
data  string.  Therefore,  it  is  necessary  to  look  at  the  upper '  lobe  peak 
amplitude  as  a  function  of  apodization  since  these  lobes  will  be  added 
to  any  real  transform  data  occurring  at  the  same  frequency.  Figure  3-4 
plots  the  relative  upper  lobe  peak  amplitude  as  a  function  of  M.  Taking 
all  three  factors  into  account,  a  bounded  optimum  value  of  M  =  30  is 
obtained,  based  on  an  upper  limit  of  2Z  peak  error  and  5Z  dc  noise. 


FIGURE  3-3. 
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4.  RESULTS  OF  COHFABATIVE  ANALYSIS 


The  cloud  images  of  che  cases  chosen  from  the  GBM  data  base  were 
analyzed  to  determine  the  results  given  in  this  report.  The  example 
scanned  f ieids-of-view  shown  in  che  illustrations  and  described  below 
represent  a  set  of  the  extremes  observed  in  these  analyses.  These  data 
were  studied  in  the  image  plane  using  the  cross-correlation  function  and 
in  the  Fourier  plane  using  power  spectral  densities  and  cross-spectra 
functions. 

4.1  ORIGINAL  CLOUD  SCENE  DATA 

The  original  cloud  data,  once  digitized,  can  be  manipulated  by 
the  computer  and  displayed  in  many  ways.  Typically,  the  GBM  mission 
data  has  been  viewed  in  Che  form  of  a  three-dimensional  plot  in  which 
the  voltage  output  from  each  detector  is  displayed  as  a  section  of  a 
plane.  These  planes  are  then  plotted  and  overlayed  into  a  three- 

dimensional  representation  of  the  scanned  field  of  view.  Figures  4-1 
and  4-2  are  three-dimensional  plots  of  Che  clouds  scenes  chosen  as 
examples  for  this  report.  Another  mode  of  display  more  nearly 
applicable  for  che  purposes  of  illustrating  clutter  data  is  che  contour 
plot.  Figures  4-3  and  4-4  illustrate  the  same  cloud  data  in  a  seven- 
level  contour  plot.  Additionally,  the  computer  allows  che  plotting  of 
the  detector  output  level  along  any  line  passing  through  the  field  of 
view.  Figures  4-5  and  4-6  are  such  plots  of  longitudinal  cuts  through 
che  contour  plots,  as  indicated  by  che  labeled  lines.  This  technique 
greatly  reduces  the  time  and  effort  required  to  analyze  these  clutter 
scenes. 


4.2  CROSS-CORRELATIONS 

The  nine  cases  containing  cloud  scenes  were  compared  between  the 
LWIR  and  SWIR  bands  in  the  image  plane  through  the  cross-correlation 
functions,  even  though  clouds  were  observed  only  in  the  LWIR  passband  in 
all  but  Case  60.  The  correlations  were  performed  between  zero  lag  and  a 
positive  lag  of  50.  This  was  done  so  Chat  any  misalignment  between  the 
two  passbands  could  be  detected  by  a  shifted  correlation  peak.  The 
solid  line  in  Figures  4-7  through  4-10  represents  the  cross-correlation 


FUiURt  /1-1.  3-d  plot,  case  49,  LWIR,  SFOV  603 
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FIGURE  4-3.  CONTOUR  PLOT.  CASE  49,  LWIR,  SFOV  603 


CONTOUR  RIOT.  CASE  60.  IWIR.  SFOV  131 
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FIGURE  4-7.  (U)  CROSS-CORRELATION.  CASE  49,  SFOV  603 
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FIGURE  4-8.  CROSS-CORRELATION,  CASE  60,  SFOV  131 


between  the  two  passbands  for  a  single  detector  output.  The  dashed  line 
illustrates  the  same  cross-correlation  except  that  the  output  of  three 
adjacent  detectors  in  the  SFOV  have  been  averaged.  The  dash/dot  line 
and  the  dash/double-dot  line  represent  cross-correlations  for  five  and 
seven  detectors,  respectively.  The  purpose  of  this  averaging  is  to 
determine  the  effect  on  the  correlation  level  of  a  varying  detector 
size,  or  a  reduced  optical  resolution.  It  is  noted  that  for  the 
uncorrelated  cases,  the  combination  of  additional  detectors  produced 
significant  variations  in  the  cross-correlation  function,  while  for  Case 
60  little  change  was  seen.  This  is  accounted  for  by  the  fact  that  the 
low  spatial  frequencies  of  the  cloud  scenes  are  the  major  contributors 
to  the  correlation  level  in  the  majority  of  cases.  This  effect  will  be 
discussed  in  more  detail  in  Section  4.3. 

Case  60  produced  significant  levels  of  correlation  between  bands 
similar  to  the  example  given  in  Figure  4-8  in  approximately  87%  of  the 
frames.  Figure  4-9  illustrates  an  interesting  negative  correlation 
which  was  observed  in  one  portion  of  the  sky  of  Case  60.  Figure  4-10 
illustrates  a  phase  lag  of  approximately  40  sampled  units  exists  between 
the  LWIR  and  SWIR  radiometric  profiles.  It  has  been  determined  that 
this  is  not  due  to  a  misalignment  of  the  two  SFOVs  as  was  first  sus¬ 
pected,  but  rather  is  a  phenomenon  of  the  cloud  itself.  These  scenes 
and  possible  scenes  showing  negative  correlation  are  indicative  of  a 
shift  in  radiometric  profile  in  clouds  between  the  LWIR  and  SWIR  pass- 
bands. 

4.3  POWER  SPECTRAL  DENSITIES  AND  CROSS-SPECTRA 

The  power  spectral  densities  of  single-  and  multiple-detector 
combinations,  as  shown  in  the  example  of  Figure  4-11,  were  calculated 
for  each  of  the  chosen  cases;  in  general,  the  frames  in  which  clouds  are 
present  exhibit  1/f  type  spectral  densities  in  agreement  with  previous 
investigations.  That  is,  the  relative  amplitude  of  the  power  spectral 
density  curve  is  inversely  proportional  to  the  spatial  (or  temporal) 
frequency. 

This  spectral  analysis  has  shown  that  in  those  cases  for  which 
olouds  appear  in  both  passbands,  the  LWIR  is  of  a  higher  overall 
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amplicude  while  Che  3WIR  axhibirs  a  greater  high-spatial-frequencv 
content.  In  chose  cases  for  which  distinct  clouds  are  visible  onlv  in 
the  LWIR,  the  noise  level  in  the  SWIR  band  in  the  region  on  the  cloud 
did  increase  about  a  factor  of  3  above  the  clear  ambient.  As  the 
density  of  Che  cloud  increases,  the  SWIR  begins  to  perceive  many  point- 
source  spikes  that  are  essentially  uncorrelated  in  time  to  the 
observable  limit  for  the  GBM  data  (i.e.,  15  Hz).  At  still  greater  cloud 
densities,  isolated  spots  of  the  cloud  formation  begin  to  grow  in  the 
SWIR  band  until  a  total  cloud  is  detectable.  Unfortunately,  this 
phenomenon  is  reduced  in  most  of  the  GBM  mission  due  to  severe 
attenuation  in  the  SWIR  passband. 


The  cross-spectra  analyses  have  indicated  that  the  major  portion 
of  the  correlation  occurs  at  the  lower  spatial  frequencies.  As  expect¬ 
ed,  the  overall  shape  of  the  cloud,  when  observed  in  both  bands,  is 
similar,  but  fine  structure  of  the  clouds  is  totally  independent  between 
bands.  This  is  a  most  important  result  since  it  suggests  that  only 
empirical  models  of  atmospheric  clutter  can  be  made  to  adequately  repre¬ 
sent  the  same  cloud  in  both  the  SWIR  and  the  LWIR  passbands,  and  that  no 
deterministic  models  can  be  developed  which  will  predict  the  radiometric 
profile  of  a  cloud  in  one  passband  when  given  the  radiometric  profile  of 
the  cloud  in  the  other. 
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i.4  ESTIMATION  OF  CLOUD  TRANSMISSION  PROFILES 

The  GBM  mission  was  to  record  the  radiometric  signatures  of 
reentry  vehicles  as  they  passed  through  the  atmosphere.  A  detailed 
model  exists  which  was  used  to  determine  the  expected  radiometric  signa¬ 
ture  of  a  heated  RV.  In  addition,  actual  measurements  of  identical  RV 
signatures  from  clear-sky  missions  can  be  accessed.  All  this  leads  to 
the  use  potential  of  the  reentry  vehicle  as  a  calibrated  IR  source.  The 
missions  for  which  the  RV  passes  behind  clouds  (there  are  many)  provide 
a  unique  opportunity  to  measure  the  transmission  profiles  of  example 
clouds.  Figures  4-12  and  4-13  illustrate  this  capability.  The  solid 
lines  represent  the  RV  irradiance  in  each  band  as  measured  by  the  GBM 
sensor,  while  the  open  boxes  trace  the  measured  RV  signature  in  clear- 
sky  conditions.  As  in  the  cross-spectra,  the  SWIR  shows  a  reduced  over¬ 


all  degradation  but  a  higher  spatial  frequency  variation  in  signal. 
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Similar  measurements  could  be  used  to  develop  a  cloud  attenuation  data 
base  from  the  GBM  mission  archives. 
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5.  CLOUD  SCENE  SYNTHESIS 


5.1  OiACT  MODEL 

The  generation  of  synthetic  imagery  to  simulate  measured  infra¬ 
red  cloud  imagery  requires  the  use  of  a  suitable  mathematical  model. 
Since  cloud  structure  is  highly  variable  and  unpredictable,  it  is  natu¬ 
ral  to  choose  for  this  purpose  some  type  of  stochastic  model.  Although 
clouds  are  inherently  three-dimensional  objects  and  three-dimensional 
stochastic  models  might  be  considered  for  their  representation,  it  seems 
more  efficient  to  attempt  to  represent  the  two-dimensional  cloud  images 
directly  rather  than  to  obtain  them  by  projection  of  three-dimensional 
cloud  structure  onto  a  plane.  This  approach  affords  a  choice  from  among 
a  large  number  of  two-dimensional  stochastic  models  which  have  been 
developed  for  purposes  of  texture  analysis  and  texture  synthesis  in  the 
fields  of  digital  image  processing  and  computer  image  generation. 

The  literature  of  texture  analysis  and  texture  synthesis  was 
surveyed  at  some  length  in  an  effort  to  identify  a  stochastic  model 
which  combines  the  attributes  of  simplicity  of  implementation  and  ade¬ 
quate  fidelity  of  modeling.  Since  these  requirements  are  obviously 
mutually  antagonistic,  a  compromise  is  to  be  expected  in  selection  of 
the  model. 

Study  of  the  literature  suggested  that  a  simple  two-dimensional 
autoregressive  model  should  satisfy  the  requirements,  particularly  since 
it  has  been  used  with  success  by  others  to  model  cloud  imagery  (Ref.  P. 
For  this  model  (as  for  most  others),  the  image  is  regarded  as  a  two- 
dimensional  rectangular  array  of  discrete  elements  (pixels),  each  of 
which  is  described  by  a  pair  of  Indices  specifying  its  location  within 
the  array  and  a  value  corresponding  to  the  image  "brightness"  or  "gray 
level"  at  the  point,  represented.  Using  the  notation  Uj^^j  to  represent 
the  gray  level  of  the  pixel  in  the  itb  row  and  jth  column  of  the  array, 
the  general  first-order  autoregressive  model  may  be  written  as  ' 

“i,j  •  *l'*i-l,j‘*'«2Ui,j-l*«3«i-l, j-1'*’  «i,j  (5-1) 


5-1 
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wh-?r- 


wh-?r-  a  . 12  *3  constants,  and  is  a  white  noise  c.  j-is-  such 

that  ^  ^  ^  ^  ^  -  •  Here  J  is  a  constant,  and  6„  ^ 

represents  the  Kronecker  delta  function. 

The  constants  ai,a2  and  aj  determine  the  autocorrelation  charac¬ 
teristics  of  the  process  defined  by  ui^j.  In  this  study,  it  is  assumed 
that  the  structure  of  cloud  imagery  is  approximately  isotropic.  Al¬ 
though  the  simple  autoregressive  model  of  Equation  1  is  inherently  ani¬ 
sotropic,  it  can  be  made  to  approximate  the  isotropic  ideal  by  choosing 
aj  *  a2  *  p  with  33  *  “  P^.  In  this  case,  0  *  1  -  p^.  Thus,  a  quasi¬ 
isotropic  form  of  the  model  of  Equation  5-1  is 

“i,o“  '“'i-l.f  'l.j  (5-2) 

In  this  model,  and  ^i,j  exhibit  the  same  variance. 

It  can  be  shown  that  the  process  generated  by  the  model  of 

Equation  5-2  has  an  autocovariance  function  of  the  form 

‘  -k,j*  (5.3) 

wher-.  denotes  the  expected  value  of  the  enclosed  random  variable, 

and  P  =E.‘-'^  ,  independent  of  i  and  j.  The  corresponding  spectral 

dens i tv  function  is 


;  ^  ,  T},  = 
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5.2  IMAGE  ANALYSIS 

To  employ  the  model  of  Equation  2  to  synthesize  imagery  repre¬ 
sentative  of  a  given  cloud  type,  it  is  necessary  to  supply  numerical 
values  for  the  model  parameters  p  and  0  defined  in  the  previous  sec¬ 
tion.  These  values  are  most  easily  determined  in  the  present  instance 
by  empirical  means  using  the  available  measured  data.  It  is  necessary 
only  to  compute  values  of  the  autocovariance  function  R(k,0)  as  defined 
in  Equation  3  for  a  range  of  values  of  the  lag  parameter  k  and  then  to 
select  values  for  p  and  o  which  give  the  best  fit  of  Equation  3  to  the 


computed  autocovariance  function.  The  value  of  o2  be  read  off 

directly  as  r(O.O).  and  p  may  then  be  identified  with  R(1,0)/  Al¬ 

ternatively,  the  value  of  p  may  be  taken  to  be  exp  (-1/k'),  where  k'  is 
the  value  of  the  lag  such  that  R(k',0)/a2  »  exp(-l). 

The  autocovariance  function  R(k,0)  was  computed  on  selected  GBM 
images,  suitably  filtered  to  resemble  data  as  it  would  appear  if  record¬ 
ed  by  the  IRST  sensor.  (See  Section  _ ).  For  this  purpose  a  stand¬ 

ard  statistical  subroutine  (FTAUTO)  from  the  IMSL  Library  of  mathemati¬ 
cal  applications  routines  was  employed.  The  autocovariance  of  the  sig¬ 
nal  from  each  of  several  selected  detectors  was  computed  for  each  detec¬ 
tor  separately,  and  then  an  average  over  detectors  was  performed  to 
yield  a  composite  value  for  R(k,0),  k  ■  1,2,..., 80.  The  resulting  com¬ 
posite  values  were  plotted  against  the  lag,  k.  The  plot  was  then  used 
to  read  off  values  of  p  and  a  as  outlined  above. 

5,3  IMAGE  SYNTHESIS 

Inspection  of  the  autoregressive  model  of  Equation  5-2  reveals 
that  samples  of  Che  process  ui^j  may  be  generated  recursively  provided 
that  values  of  j  *  1,2,...J  and  u^^i,  i  *  1,2,... ,I  are  available. 

That  is,  given  values  of  uj^^j  for  all  pixels  in  the  first  row  and  the 
first  column  of  the  image  array,  the  values  of  all  remaining  pixels  may 
be  directly  determined  by  successive  application  of  Equation  5-2.  The 
remaining  problem  is  Chen  to  provide  values  for  pixels  in  the  first  row 
and  first  column  of  Che  array. 

The  approach  taken  in  this  study  is  to  supply  the  required  pixel 
values  from  a  combination  of  measured  data  and  synthesized  values  as 
follows.  A  typical  sample  of  data  from  a  single  GBM  detector  (suitably 
filtered  .to  simulate  IRST  data)  consisting  of  a  sequence  of  93  succes¬ 
sive  pixels  is  used  as  Che  input  "time  series"  to  SUBROUTINE  FTCMP  of 
the  l^L  Library.  This  routine  performs  an  analysis  of  a  time  series 
based  on  a  stochastic  autoregressive  integrated  moving  average  (ARIMA) 
model.  It  develops  values  of  required  model  parameters  from  the  given 
tinw  series  and  then  uses  these  values  to  synthesize  a  continuation  of 
the  stochastic  process  for  future  time.  With  the  93-pixel  sample  as 
input,  FTCMP  was  used  to  extrapolate  the  sequence  to  give  a  total  of  512 


samples.  These  512  samples  were  then  used  as  the  values  of 
=•  i,2,...,512  and  again  as  values  of  u^^i,  i  =  1.2,. ..,512.  A 
straightforward  application  of  Equation  2  using  values  of  p  and  a  Dotainec 
as  described  in  the  preceding  section  then  permitted  filling  of  the 
remainder  of  the  array. 

The  procedure  described  above  may  be  used  to  synthesize  homo¬ 
geneous  cloud  images,  but  of  greater  interest  in  the  present  application 
are  celestial  scenes  containing  a  mixture  of  cloud  and  clear  sky.  The 
approach  used  in  this  study  to  generate  such  images  is  sometimes 
referred  to  as  the  Background-Foreground  model.  In  this  model,  the 
image  array  is  suitably  partitioned  into  regions  of  two  types.  Although 
Che  terms  "background"  and  "foreground"  are  conventionally  used  Co  label 
these  regions,  "sky"  and  "cloud"  will  serve  better  here.  Homogeneous 
images  of  sky  and  cloud  are  generated  as  described  above,  and  the  par¬ 
titioned  array  is  used  as  a  guide  in  constructing  a  composite  image 
using  Che  following  simple  algorithm.  Pixels  in  the  partitioned  image 
are  examined  sequentially,  and  their  type  ("sky"  or  "cloud")  is  noted. 
If  the  pixel  type  is  "sky,"  the  value  of  the  pixel  at  the  corresponding 
location  in  Che  homogeneous  sky  image  is  inserted  into  the  composite 
image  at  the  corresponding  location.  Similarly,  if  the  type  "cloud"  is 
noted,  the  pixel  at  the  corresponding  location  in  the  composite  image  is 
assigned  the  value  of  the  pixel  at  Che  corresponding  position  in  the 
homogeneous  cloud  image. 

Use  of  the  Background  Foreground  approach  introduces  two  new 
problems:  How  is  the  partitioning  to  be  accomplished  to  give  a  realis¬ 

tic  (irregular)  boundary  line  separating  "cloud"  and  "sky"  regions?  How 
can  Che  abrupt  transition  which  will  occur  at  the  boundary  be  made  more 
realistic  without  additional  processing? 

The  partitioning  problem  is  solved  by  employing  stochastic 

interpolation  schemes  proposed  by  Fournier,  Fussell,  and  Carpenter 
(Ref.  2)  to  approximate  fractal  curves.  It  has  been  shown  (Ref.  3)  that 
cloud  perisieters  are  fractal  curves  of  fractal  dimension  0  >  1.35.  The 
method  of  Reference  2  permits  the  generation  of  a  fractal  curve  of  the 

appropriate  fractal  dimension  between  any  pair  of  points  in  a  plane.  In 

the  present  application,  one  of  these  points  is  chosen  Co  lie  on  the 
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upper  boundary  of  the  image  plane,  while  the  second  point  is  chosen  to 
lie  on  the  lower  boundary.  The  fractal  curve  connecting  these  two 
points  thus  divides  the  image  roughly  into  left  and  right  regions  which 
may  be  identified  with  "sky"  and  "cloud",  respectively.  A  slight  modi¬ 
fication  of  the  method  of  Reference  2  permits  the  introduction  of  a  con¬ 
trolled  degree  of  convexity  into  the  generated  curve  to  simulate  the 
"fluffy"  appearance  of  clouds.  Without  this  feature,  anv  curves  gene¬ 
rated  would  be  statistically  symmetrical,  and  the  assignment  of  "sky" 
and  "cloud"  values  to  left  and  right  regions  would  be  completely  arbi¬ 
trary.  The  convexity  parameter  introduced  is  simply  a  constant  bias 
value  added  to  the  random  variable  used  in  the  stochastic  interpolation 
algorithm  of  Reference  2.  Since  no  empirical  value  for  this  bias  has 
been  reported,  a  value  (0.02)  was  selected  which  appears  visually  to 
yield  credible  cloud  contours. 

Plots  of  infrared  cloud  data  for  scans  which  intersect  cloud 
boundaries  typically  exhibit  a  gradual  decrease  in  cloud  radiance  as  the 
cloud-sky  boundary  is  approached  from  the  "cloud"  side.  In  addition,  as 
the  cloud  becomes  more  tenuous  near  the  boundary,  there  may  be  "holes" 
in  the  cloud  through  which  the  sky  may  be  seen.  To  represent  these 
effects,  the  "cloud"  portion  of  the  image  is  modified  as  follows. 

Let  C(x,y)  and  S(x,y)  represent  the  cloud  and  sky  radiances  at 
point  (x,y)  in  the  two  homogeneous  regions,  respectively.  Consider  a 
single  horizontal  scan  line  (row  of  pixels)  from  the  image  corresponding 
to  a  constant  value  of  y,  and  suppose  that  the  cloud  boundarv  intersects 
this  scan  line  at  the  point  (xj,y)  with  "cloud"  corresponding  to  x^  , 
"sky"  to  X  <xi.  We  define  a  modified  cloud  radiance  C'(x,y)  for  x  >  x: 
by  C'(x,y)  =•  greater  of  (S(x,y)  and  C  (x,y)  -  C(xi,y)  exp  ( -(x-xj ) /d) . 
Here,  d  is  a  characteristic  length,  determined  empirically,  which  pro 
vides  a  measure  of  the  rate  of  decrease  of  cloud  radiance  as  the  cloud- 
sky  boundary  is  approached.  In  the  case  where  the  horizontal  scan  line, 
y,  intersects  the  boundary  at  two  points,  xj  and  X2,  with  the  interval 
(XpX2)  corresponding  to  "cloud",  the  expression  for  C'(x,y)  with  x^  ^  x  < 

X2  is 

C'(x,y)  - 

greater  of  (S(x,y)  and  C(x,y)-C(xj  ,y)exp(-(x-x^) /d)-C(x2  ,y)exp(-(x.,-x) /d) )  . 


Figures  5-1  through  5-7  are  graphic  representations  of  cloud 
scenes  generated  using  the  above  method.  Figures  5-1  and  5-2  represent 
clouds  generated  based  on  the  data  for  Case  60  for  the  LWIR  and  SWIR 
passbands,  respectively.  Figure  5-3  is  also  based  on  Case  60  SWIR  data 
but  for  a  reduced  sensor  sensitivity  and  resolution. 

5.4  TEMPORAL  PROPAGATION 

In  addition  to  the  ability  to  synthesize  textures  resembling  the 
spatial  distribution  of  infrared  radiation  observed  in  clouds,  it  is 
also  desirable  to  be  able  to  simulate  the  evolution  of  these  distribu¬ 
tions  in  time.  Given  a  synthetic  radiance  distribution  at  time  t  =  ta 
represented  by  an  array  of  pixels  a^^j,  such  a  model  would  be  capable  of 
providing  a  credible  distribution  b^^j  over  the  same  set  of  pixels  at  a 

later  time  t  =  tj,.  Any  model  proposed  for  generating  the  array  b£  j 
should  exhibit  certain  statistical  characteristics.  Thus,  in  general, 
the  array  b^^j,  while  differing  in  a  random  manner  from  a^^j,  should 
nevertheless  be  correlated  with  ai^j  to  a  degree  dependent  on  the  value 
of  the  time  difference  t^  -  t^.  In  particular,  when  t  =  tg,  the  model 
must  predict  =  ai^j  for  all  i  and  j.  For  t^  greater  than  tg,  the 

model  should  predict  a  monotonically  decreasing  value  of  the  quantitv 
Hfai^jbi^j;  as  the  difference  t>j  -  tg  increases,  and  in  the  limit  as  tij 
approaches  infinity,  the  expected  value  E  [ai^jbi^j]  should  approach 
zero.  In  addition,  the  spatial  covariance  properties  of  the  b-arrav 
should  be  identical  to  the  corresponding  properties  of  the  a-array  for 
all  values  of  tj,  and  tg,  at  least  for  time  differences  t^j  ~  tg  on  the 
order  of  a  second.  In  particular,  we  must  have  E  [  b^^  j  ]  =  E  [  a^  j  ]  and 
E  t*’i,  jt>i+h,  j+k  1”  ^  j^i+h,  j+k^  integral  values  of  i,  j,  h,  and 

k. 

Consider  a  model  for  b^^j  of  the  form 


V  \  X  ’.  •V'.  ' 


The  parameter  p  is  a  function  of  the  time  t^at  which  the  distri¬ 
bution  is  given  by  The  array  may  be  thought  of  as  the  radi¬ 

ance  distribution  existing  after  an  infinite  time  has  elapsed  following 
observation  of  the  distribution  at  time  tg.  A  reasonable  model  for 

the  dependence  of  p  on  the  time  difference  t^j  -  tg  is 

?  ('^b  "  ^a)  “  exp  L  -(tb  -  Ca)  '*-'. 

where  "  is  the  characteristic  time  associated  with  the  evolution  of  the 
radiance  distribution.  The  value  of  r  could  be  determined  empirically. 
With  this  expression  for  p,  the  model  for  becomes 

bi^j  -  exp  (-(tb  -  ta)/T)  aij 

♦  ^1  -  exp  (-2(tb  ■  La^/r)  ci^j 

which  satisfies  the  requirements  discussed  above. 

where  E  [  ai^j]  *  0,  and  c^^j  is  a  new  array  generated  using  the  same 
stochastic  process  as  that  used  in  generating  the  a-array  but  statis¬ 
tically  completely  independent  of  the  a-array.  That  is, 

j^i+h,  j+k- 


for  any  i,  j,  h,  and  k,  but 


f  ®i,  j‘^i+h,  j+k  1  ® 


for  all  i,  j,  h,  and  k.  Then  we  have 


®  Ibi^ jbi+h, j+k^  *  t*i, j*i+h, j+kl 

+  (l-p2)  E  [ jci+b, j+k  1  "  ®  [ *i, j^i+h, j+k)  • 
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APFERDIX  A.  DATA  lEDncnOIl 


This  appendix  outlines  the  salient  features  of  the  data  reduc¬ 
tion  procedures  only  to  the  extent  required  to  appreciate  the  signifi¬ 
cance  of  the  data  presented  in  the  main  body  of  this  report. 

A. 1  TIME  CONVERSION 

Throughout  the  data  reduction  process,  TALO  was  used  exclu¬ 
sively,  in  particular  to  merge  data  from  multiple  independent  data 
streams.  TALO  was  obtained  by  subtracting  launch  time  from  Greenwich 
Mean  Time  (GMT)  recorded  on  the  original  analog  tape.  GMT  times  record¬ 
ed  with  TRIG  B  were  used. 

A. 2  RECORDED  VOLTAGES 

The  detector  output  voltages  are  recorded  as  eight-bit  binary 
words,  each  of  which  can  assume  at  most  2®,  or  256,  distinct  values. 
The  quantitizing  associated  with  this  recording  is  designed  to  remain  at 
less  than  5*  over  the  system's  dynamic  range  by  use  of  a  piecewise 
linear  compression  scheme. 


A. 3  FILTERS 

The  decoded  PCM  data  from  each  detector  may  be  filtered  to 
eliminate  noise.  The  characteristics  of  the  filters  used  depend  on  the 
nature  of  the  noise. 

A. 4  PEAK  DETECTION 

Detection  of  pulses  arising  from  the  presence  of  point  sources 
in  the  field  of  view  is  based  on  fitting  detector-dependent  pulse  shapes 
derived  from  measurements  of  a  calibration  point  source.  The  method  of 
•least  squares  is  used  to  determine  three  free  parameters:  1)  the  pulse 
amplitude,  2)  the  local  baseline  level,  and  3)  the  slope  of  the  local 
baseline.  The  variance  of  the  fit  is  computed  and  is  used  as  a  detec¬ 
tion  criterion. 

When  a  local  maximum  in  amplitude  is  detected  that  exceeds  a 
preset  threshold  and  is  coincident  (to  within  a  specif led .tolerance) 


vith  a  minimum  in  standard  deviation.  a  oreliminary  detection  is 
declared.  Following  preliminary  detection,  additional  constraints  are 
imposed,  as  needed,  to  ensure  elimination  of  noise  or  radar-induced 
false  targets  while  retaining  the  true  target. 


A. 5  CON^/ERSION  OF  PULSE  AMPLITUDE  TO  PHYSICAL  UNITS 

Application  of  the  system  calibration  to  target  measurements  is 
based  on  the  assumptions  that  1)  most  of  the  energy  in  a  point  source 
image  is  distributed  over  a  detector  triplet  (three  adjacent  detectors) 
and  2)  the  distribution  of  energy  in  the  focal  plane  is  the  same  for 
both  target  and  calibration  source  images.  When  more  than  three  adja¬ 
cent  detectors  record  pulses  from  a  single  target,  then  the  detector 
triplet  selected  is  the  one  for  which  the  combined  pulse  is  the  strong¬ 
est.  However,  if  fewer  than  three  detectors  record  pulses  (adjacent 
detector(s)  having  zero  amplitude),  the  detector  registering  the  larger 
pulse  is  selected  as  the  center  detector  of  the  triplet. 

The  three  amplitudes  of  the  target  pulse  are  compared  to  each 
triplet  of  the  same  detector  set  in  the  calibration  table  (refer  to 
Appendix  G) ,  and  a  calibration  triplet  is  selected.  The  selection  is 
based  on  how  closely  the  amplitudes  of  the  calibration  triplet  resemble 
the  target  triplet  of  the  same  detector  set.  This  criterion  ensures 
that  the  calibration  applied  to  the  target  measurements  is  the  one  for 
which  the  distribution  of  calibration  irradiance  over  the  three  detec¬ 
tors  most  closely  matches,  within  measurement  uncertainties,  the  distri¬ 
bution  of  target  irradiance  over  the  same  set  of  detectors. 

The  effective  irradiance  of  the  target  at  the  GBM  system  is 
estimated  by  the  ratio  of  the  measured  target  triplet  pulse  to  the 
selected  calibration  triplet  pulse  times  the  calibration  source  irradi¬ 
ance  at  the  time  of  the  recorded  calibration. 

% 

A. 6  COMPUTATION  OF  RADIOMETRIC  PARAMETERS 

A. 6.1  Intensity  from  Irradiance 

The  spectral  irradiance  at  an  aperture  by  a  source  can  be  given 
(X.  T)  -  N,v^(^,  T)  T(X)  .  (A-1) 


where 


a  ~  projected  area  of  the  source,  assumed  to  be 

uniformly  heated  to  temperature  T 

R  ■  line-of-sight  range  to  the  source  (dimensions 

of  length) 

£(a)  ■  emissivity  of  the  surface  of  the  source, 

assumed  uniform  over  the  area  a 
(dimensionless) 


”  spectral  radiance  of  a  blackbody  at  wavelength 
X  and  temperature  T  as  given  by  the  familiar 
theoretical  expression  of  Planck  (dimensions 
of  power  •  area*'  •  solid  angle~^  •  wavelength”^) 

t(X)  «  transmission  along  the  line  of  sight  as  a 

function  of  wavelength  (dimensionless) 

(U)  The  effective  in-band  irradiance  (in  either  band,  b).j^ 

^eff  ’ 

is  defined  as  ° 


y*  T)  'r^(X)  dX  ,  (A-2) 

0 

where  is  a  system  normalized  responsivity  of  band  b  at  wavelength 

X.  (This  is  the  total  spectral  responsivity  of  the  detector  as  installed 
in  the  system  and  includes  the  effect  of  the  detector,  optics,  filter, 
etc.)  The  in-band  radiant  intensity  from  a  source  is  defined  as 


*R^(X)  dX 


^  A“3 


where  J, (X,  T)  is  the  spectral  radiant  intensity  of  the  source  at  wave- 

A 

length  X;  and  the  spectral  irradiance  is 


T(X) 

'P 


(A-4) 


Therefore,  by  cottblnlng  Equations  A-3  and  A-4, 


Since  GBM  is  a  two-color  system  and  cannot  make  spectral 
measurements,  the  effective  in-band  irradiance,  is  extracted  using 

the  method  described  In  Section  A. 5  rather  than  being  measured  directlv. 

A 

In  order  to  evaluate  Equation  A-3 ,  t(X),  as  weighted  by  and  must 

be  estimated.  The  atmospheric  transmittance  t.  (3)  is  estimated  bv  a 

—  D 

wavelength- independent  quantity  *  (T)  defined  as 

b 


where 


jT  (X.  T)  T(X)  ‘R^(X)  dX 

7  (T)  =  ■  ■  -g -  • 

b 

J  (X.  T)  (X)  dX 

0  p 

J,  (X,  T)  ■  spectral  radiant  Intensity  of  the  source  as 


uA-6) 


p 

T(X) 


predicted  by  the  OSC 

atmospheric  transmission  as  calculated  using 
LOWTRAN  4 


By  substituting  T  for  t(X)  in  Equation  A-5 
J  m  ■  ^ \(X)  d-.  . 


and  by  substituting  Equation  A-S,  the  working  equation  is  obtained: 
^ef  i. 


A--) 


^  ... 
err. 


A . b . 2  Two-Color  Parameters 

Calculation  of  the  the  two-color  parameters,  which  are  two- 
color  temperature  (T) ,  emissivity  area  (ea) ,  and  broadband  intensity  (J) , 
is  described  below. 

The  two-color  temperature  is  calculated  from  the  ratio 
■'.ff  /  T)  ■' 


r(T)  - 


B  0 


‘Rj(X)  dX 


•^eff  /** 

J  J(X.  T)  ‘ 
0 


(A-9) 


R^(X)  dX 


A-4 


I 


wh«r«  Che  subscript  B  refers  to  the  SVIR  band  and  Che  subscript  A  refers 
CO  the  LWIR  band.  defined  in  Zauation  A-3.  Substituting 

T)  -  a  T)  in  Equation  A-3  gives 


i 

I 

i 

s 

i 

e 


j 


e£f 


T) 


^R.  (X)  dA 
0 


(A- 10) 


and  Equation  C-9  becomes 


r(T) 


■A 


(X)  T) 


(X)  N^jjCX,  T) 


dX 


^R»(X)  dX 

A 


(A-ll) 


There  are  four  unkno%nis  in  Equation  A-ll:  the  emitting  areas  of 
the  target  in  the  LWIR  and  SWIR  bands,  a^,  and  ag,  and  the  spectral  emit- 
tances  in  each  band.  If  it  is  assumed  that  the  emittance  (\)  is  wave¬ 
length  independent  and  that  a^  =  ag,  and  the  appropriate  cancellations 
are  made  in  Equation  A-ll, 


r(T) 


T)  1r„(X)  dX 


/  \b('. 


T)  ^R^(X)  dX 


(A-12) 


The  integrals  appearing  in  Equation  A-ll  were  evaluated  numerically 
using  16-point  Gaussian  quadrature.  Figure  A-1  is  a  plot  of  r(T)  as  a 
function  of  temperature. 

One  determines  the  color-temperature  from  the  mission  data  by 
associating  the  effective  radiant  intensities  (obtained  from  the  meas¬ 
ured  irradiances  by  means  of  Equation  A-8)  with  the  temperature  in 
Figure  A-1. 


J: 


SS 


The  emissivity  area  (a  €  )  is  found  using  the  two-color  tempera¬ 
ture  and  Equation  A-10.  One  further  assumes  that  the  emissivity  C'X)  is 
independent  of  the  wavelength.  Solving  Equation  A-10  for  the  emissiv¬ 
ity  area  gives 


A-5 


SO  □ 

where  either  waveband  may  be  used,  the  results  being  identical  in  the 
two  cases. 

The  total  broadband  radiant  intensity  of  the  source  is  computed 
using  the  Ste f fan-3ol tzmann  law  expressed  by 

»  (A-14) 


J  - 


ae  a  T 


IT 


where  o  is  the  Steffan-Boltzmann  constant,  aC  is  given  by  Equation  A-13. 
and  T  is  the  two-color  temperature. 

A. 6. 3  Smoothing  Techniques 

The  effective  irradiance  and  the  effective  radiant  intensity  are 
smoothed  by  a  running  average  taken  over  a  0.5-sec  interval  of  data  (15 
adjacent  SFOVs).  For  this  average  to  be  calculated,  at  least  eight 
SFOVs  must  contain  valid  data.  For  each  interval,  the  averaged  data  are 
assigned  the  altitude  of  the  midpoint  of  the  interval.  This  interval  is 
then  stepped  one  SFOV  at  a  time,  and  the  procedure  repeated. 

A. 6. 4  Standard  Deviation 

For  single-color  parameters  and  J  ),  the  standard 


deviations  are  defined  as 


eff,. 


N  -  1 


efl 


N  -  1 


(A- 15 


(A-16') 


where 


*  mean  value  of  effective  Irradiance 
b 

*  mean  value  of  effective  radiant  intensity 


N 


number  of  points  being  averaged. 


For  two-color  parameters,  the  standard  deviations  are  defined  as 


-  f  +  Q(T)  (A-17) 

where  f  is  the  function  referred  to  in  Equation  A- 10  and 


Q 


(A-18) 


This  method  of  defining  a-£  propagates  the  oJeff^  errors  under 
a  worst-case  assumption.  Similarly,  worst-case  error  propagation  is 
used  in  calculating  the  errors  of  Og  and  Oj. 


A. 7  PRECISION  AND  REPEATABILITY 


Definitions  of  precision  and  repeatability  have  been  adapted  to 
GBM  data  as  follows:  precision  denotes  the  degree  of  scatter  in  each 
mission's  data  and  is  defined  at  a  given  altitude  as 


where 

Xi 

u 


-  1 


(A-19) 


=  the  measured  radiometric  values 

=  the  average  value  of  over  the  0.5-sec  smoothing 
interval 


N 


»  the  number  of  measurements  during  the  smoothing 
interval . 


Precision  is  expressed  in  percent  and  is  the  average  of  the  precision  at 
each  altitude  for  the  entire  LWIR  or  SWIR  data  span: 


A-1  " 


(A-20) 


where  R  is  the  number  of  altitudes  for  which  is  calculated. 


Repeatability  denotes  the  degree  of  disagreement  between  data 
sets  (missions)  separated  by  a  significant  amount  of  time  and  is  defined 


at  a  given  altitude  as 


where 


^  -  1** 
X  (u  -  u)^ 
.1»1  ^ 

M  -  1 


(A-21) 


=  the  smoother  radiometric  data  of  mission  i  at  a 
given  altitude 

=  the  average  of  ui  over  all  missions 
=  the  number  of  missions  (M^2). 


Repeatability  Is  expressed  In  percent  and  is  the  average  of  the  repeat¬ 
ability  at  overlapping  altitudes: 


L  P. 

P  •  "Z  -  ^00/1 

A-1  u 


(A-22) 


where  L  is  the  number  of  altitudes  for  which  overlapping  data  (data  taken 
at  altitudes  of  mutual  coverage)  exist  (L  ^  2) .  These  altitudes  of  mutual 
coverage  have  been  interpolated  to  even  altitudes  at  2-km  intervals.  For 
data  gaps  greater  than  twice  the  smoothing  interval,  the  interpolated 
value  of  p.  has  oeen  omitted. 


APPENDIX  B.  CALIBRATION 


3.1  MEASUREMENTS 

Following  the  mission,  and  prior  to  securing  the  GBM  Sensor  Svs- 
tem  from  its  operational  cooled  status,  a  postmission  radiometric  cali¬ 
bration  is  recorded.  This  calibration  is  used  as  a  record  of  the 
detector  responsivility  for  the  mission,  and  is  performed  as  closely  as 
practical  to  the  time  of,  and  under  the  same  conditions  as,  the  mission 
being  recorded.  This  calibration  is  verified  by  use  of  a  premission 
calibration. 

The  radiometric  calibration  is  performed  using  the  calibration 
collimator  designed  especially  for  this  purpose.  Basically,  the  colli¬ 
mator  consists  of  a  graybody  source  controllable  over  the  range  from 
about  300  to  1200  K,  a  set  of  apertures  that  may  be  placed  in  front  of 
the  graybody.  and  a  system  of  collimating  optics  that  present  to  the 
sensor  a  virtual  source  at  infinite  range.  While  the  GBM  tracking 
pedestal  is  held  at  the  central  azimuth  location  and  the  GBM  primarv 
telescope  is  moved  in  elevation,  the  graybody  image  is  slowly  stepped  up 
or  down  the  focal  plane  over  the  elevation  field  of  view. 

The  calibration  coefficients  have  been  derived  using  these  data 
and  various  assumptions.  Many  of  the  assumptions  have  been  verified  to 
within  limits  by  analysis  of  more  extensive  calibrations,  such  as  those 
that  were  performed  during  laboratory  tests  prior  to  shipment  of  the  GBM 
system  to  Roi  Namur  for  installation. 

One  assumption  used  in  calibration  is  that  the  GBM  system  is  a 
linear  system,  since  evidence  accumulated  thus  far  suggests  that  the 
sensor  response  is  nearly  linear  over  its  full  dynamic  range.  Another 
important  assumption  used  in  calibration  is  that  the  response  of  the 
system  to  a  point  target  can  be  described  by  a  triplet  of  voltage 
values.  This  is  apparent  from  observations  that  the  voltage  values  from 
a  point  target  are  frequently  detected  by  more  than  one  detector,  while 
hits  on  three  adjacent  detectors  are  very  common. 

B.2  THEORY 

Calibration  of  the  sensor  is  based  on  the  premise  that  one  can, 
through  a  combination  of  measurements,  calculate  the  effective 


irradiance  present  at  the  sensor.  The  calibration  method  applied 
requires  only  chat  the  energy  distributions  in  mission  and  calibration 
source  images  be  similar.  This  allows  data  obtained  with  both  focused 
and  moderately  defocused  systems  to  be  reduced.  The  method  allows  the 
total  irradiance  to  be  determined  even  when  a  portion  of  the  energv 
falls  on  an  inactive  detector  element. 


The  basic  theory  employed  in  computing  the  effective  irradiance 
produced  by  an  unresolved  source  is  based  on  the  condition  chat  the 
source  irradiance  h^,  during  a  single  calibration  run  need  only  remain 
constant  during  the  recording  of  any  voltage  triplet.  These  voltage 
triplets  thus  satisfy  an  equation  of  the  form 

''bji-  ''bU+Di^  ■  b.  •  is-i, 


where 


V.  -  measured  waveform  voltage  above  the  local 
baseline  of  the  Jth  detector  in  spectral 
band  b  of  the  1th  pulse 


h^  a  source  irradiance  at  the  sensor  in  spectral 
band  b. 


This  function has,  assuming  only  linearity  of  detector  response,  the 


property 


[5V, 


b(J-l)’ 


^b(j+l)-  "  ’^bj’  ’'^bCj+D- 


where  5  is  any  arbitrary  factor.  Choosing  5  equal  to  the  square  root  cf 
the  sum  of  the  squares  of  Che  triplet  voltage  yields 


^^^b(j-l)*  ^bj’  \a+l)^ 


where 


bCj-l)  ’bj  ^  ’b(J+l) 


®bj  “  ^b(j-l)^l^lbj 
®bj  "  ^bj^l^lbj 


^bj  ■  \(j+l)' 


m 


This  choice  of  5  permits  to  be  expressed  as  the  product  of  two  factors, 
one  of  which  depends  only  on  the  direction  cosines  in  Che  voltage  domain. 
Expressed  in  this  manner  represents  a  universal  function  that  Is  inde¬ 
pendent  of  Che  particular  Irradlance  used  in  obtaining  the  calibration 
data.  The  unknown  Irradlance  can  therefore  be  obtained  by  multiplying 


3.3  NDISE  CONSIDERATIONS 

Assuming  that  errors  in  Che  measurements  caused  by  random  fluc¬ 
tuations  (noise)  are  additive  and  normally  distributed  about  the  mean 
value,  Chen  each  member  of  a  measured  voltage  triplet  may  be  represented 
in  the  form  V  *  V  +  Av,  where  V  is  the  true  (mean)  value  and  Av  is  the 
random  error,  which  is  normally  distributed  having  the  variance  ^  • 
The  variance  may  be  different  for  each  detector  element. 

This  leads  to  the  concept  of  a  "resolution  cell"  in  the  voltage 
domain.  The  cell  is  ellipsoidal  in  shape  and  is  described  by  the  func¬ 
tion 


where  for  detector  J,  etc.  The  value  of  the  constant  in  Equation 

B-4  is  somewhat  arbitrary,  but  if  the  50%  point  is  selected,  the  value 
is  2.36606.  and  there  is  a  50%  probability  chat  the  observed  triplet 
V(j_i),  Vj,  V(j+i)  falls  within  Che  ellipsoid  so  defined. 

The  existence  of  a  resolution  cell  of  finite  size  in  voltage 
space  suggests  that,  within  the  allowed  finite  ranges  of  voltages,  there 
are  only  a  finite  number  of  distinguishable  resolution  cells.  In  par¬ 
ticular,  if  the  number  of  voltage  triplets  measured  during  calibration 
exceeds  Che  number  of  distinguishable  resolution  cells  in  voltage  space 
as  determined  by  random  noise,  it  is  justified  to  combine  some  of  Che 
measured  triplets  by  averaging  to  reduce  Che  total  number. 

Therefore,  for  each  pair  r,  s  (r  *  s,  s  »  1  to  n,  where  n  is  the 
total  number  of  voltage  triplets  for  one  detector  triplet),  the  value  of 
Q73  is  found  by 
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U-x),. 


'(j-l) 


(j+1) 


where  W  is  the  smaller  of  W  and  W  .  The  inclusion  of  the  factor  U  in 
X  r  3  X 

the  formula  for  reflects  the  shrinking  size  of  the  resolution  cell  when 
points  are  combined  and  makes  it  less  likely  that  two  points  having  large 
and  nearly  equal  weights  will  be  combined.  The  general  effect  is  to  cause 
"clusters"  of  points  in  voltage  space  to  coalesce  to  a  single  point  each, 
with  the  weights  of  the  final  set  corresponding  to  the  number  of  points  in 
each  original  cluster.  Let  denote  the  smallest  calculated  above. 

If  <  2.36606,  then  the  points  k  and  I  are  combined  by  averaging  accord¬ 
ing  to  the  formulas 


w  -  , 


’(J-l)  *  ’(J-l)k  "l  ’(J-l)l'^ 

Vj  -  (W^Vj^.WjVj^)/W  , 

’(J+1)  -  '"k  ’(J.M)k  *  «1 


(B-6) 


The  two  original  points  k  and  I  are  then  replaced  by  the  single  average 
point,  thereby  reducing  the  total  number  by  one.  The  process  is  repeated 
until  >  2.36606. 


When  no  more  points  can  be  combined,  the  resultina  triolet 
voltages  are  converted  to  direction  cosines  i.  ..  ,  which 

JJIX  OJ«t  OJ  iC 

together  with  the  voltages  constitute  the  calibration  table  for 

the  center  detector. 


B.4  CALIBRATION  FACTORS 

The  calibration  factors  were  derived  by  measuring  the  value 

of  the  waveform  voltage  above  the  local  baseline.  These  factors  were 
obtained  as  the  constant  image  of  the  collimating  graybody  source  slowly 
moved  down  each  detector  in  elevation.  Different  responsivities  were 
measured  when  the  secondary  mirror  scanned  the  image  across  the  focal 
plane  right  -to  left  and  when  it  scanned  left  to  right  (see  Figure  2-3 


for  focal  plane  detector  arrangement-.  This  direction-correlated  effect 
was  accounted  for  in  extraction  of  the  intensities. 

Normally,  mission  data  are  converted  to  engineering  units  on  the 
basis  of  calibration  data  recorded  immediately  before  and/or  after  the 
mission.  However,  in  the  case  of  mission  STM-13W.  calibration  data  were 
recorded  only  on  two  detectors  of  the  LWIR  array  and  two  detectors  of 
the  SWIR  array.  The  reduction  process  requires  a  calibration  for  each 
detector  separately;  therefore,  the  calibration  data  recorded  for  mis¬ 
sion  GT  134M  were  substituted  in  the  conversion  of  measured  voltages  to 
units  of  irradiance.  This  calibration  file  was  selected  because  it  was 
of  good  quality  and  was  close  in  time  to  STM-13W. 

To  compensate  for  the  relative  responses  of  the  detectors 
between  the  calibrations,  it  was  necessary  to  develop  and  apply  a  set  of 
correction  factors.  These  procedures  are  outlined  in  Figure  B-1.  As 
shown  in  the  figure,  the  substituted  calibration  data  were  processed  in 
the  normal  manner  to  generate  a  calibration  file.  Voltage  peaks  for  the 
mission  calibrations  were  read  by  hand  from  strip  charts  and  converted 
to  irradiances  by  means  of  the  substituted  calibration  file.  Since  the 
effective  irradiances  during  the  mission  calibrations  are  known,  the 
departure  of  the  "calibration"  irradiances  from  this  value  represents 
the  error  to  be  corrected  using  suitable  correction  factors. 

Correction  factor®  CJ  defined  by 

ci  -  H,<a;> 

were  developed,  where  i  -  1,  2,  ....  31  Identifies  the  particular  LWIR 
detector  and  ±  indicates  the  direction  of  scanning.  The  symbol 
denotes  the  average  value  of  the  ''caulibration"  Irradiance  for  the  ith 
detector  and  the  ±  scan  direction.  Finally,  the  correction  factors  are 
used  together  with  the  calibration  data  file  to  convert  the  mission 
voltages  to  irradiance  values. 


FIGURE  B-1.  '  SCHEMATIC  DIAGRAM  ILLUSTRATING  THE  METHOD  OF  DEVELOPING  AND  APPLYING  CORRECTION 

FACTORS  C!  f()  COMPENSATE  FOR  CHANGE  IN  RELATIVE  RESPONSE  OF  DETECTORS  ON  MISSION 
STM-I3W 


APPKRDIX  C.  METEOROLOGY 


C- 1  RADIOSONDE  BALLOON  FLIGHTS 

For  each  GBM  mission,  radiosonde  balloon  flights  were  conducted 
to  collect  meteorological  dat  supporting  the  Army  Optical  Station.  These 
data  included:  wind  direction,  wind  speed,  temperature,  dew  point, 
pressure,  density,  absolute  and  relative  humidity,  index  of  refraction, 
visibility,  and  wind  shear  as  a  function  of  altitude  between  12  and 
100,000  ft. 

Those  data  of  most  interest  in  the  analysis  of  clouds  are  the 
temperature  and  dew  point  profiles.  These  are  plotted  for  the  cases  of 
interest  between  0  and  40,000-ft  altitude  in  Figures  C-1  through  C-18. 


C.2  CLOUD  RANGE  ESTIMATES 

A  rough  estimate  of  the  range  to  an  observed  cloud  can  be 
determined  using  the  bores ight  tracking  data  of  Appendix  C  and  the 
Temperature  and  Dew  Point  Profiles.  These  estimates  are  determined 
assuming  that  the  temperature  and/or  dew  point  inversions  as  illustrated 
in  Figure  C-19  indicate  the  average  altitude  for  cloud  formations.  Then, 
using  Che  boresight  elevation  of  the  cloud  measurement  and  simcle 
tr icor.cmetry ,  a  rougn.  estimate  or  tne  cloud  range  is  made.  Ine  ran^e 
error  associated  with  such  an  estimate  is  aiven  bv. 


a 


(C-l 


where  the  elevation  angle  error  Og  is  insignificant  compared  to  the 
altitude  error  0^;  hence  Equation  C-l  reduces  to 
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%  ^  -i- 


(C-2) 


where 


*  altitude  error  (1  sigma) 

9  »  measured  elevation  angle  of  cloud. 


I 


C-l 


TLMPERATURt;  (°C) 

I  K.Dia  C-/.  CASE  42  TEMPERATURE  PROFIL 


vs  ALTI 


-20 

I EMPERA 
CASE  46  TEMI 


aaniiiiv 


I 


in  <9  r)  <9  in  0)  in  <9* 

ro  M  «\i  <M 


3aniuiv 


C-14 


I 


ruinjirLm:'4Vvv\.%uv  Lnt  l 


TEMPERATURE  {°C) 

FIGURE  C-13.  CASE  49  TEMPERATURE  PROFILE 


EW  POINT  (“C) 

iiE  49  DEW  POINT  vs  ALT 


RANGE  ESTIMATION 


APPENDIX  B.  CALIBRATION 


3.1  MEASUREMENTS 

Following  the  mission,  and  prior  to  securing  Che  GBM  Sensor  Svs- 
cem  from  its  operational  cooled  status,  a  postmission  radiometric  cali¬ 
bration  is  recorded.  This  calibration  is  used  as  a  record  of  the 
detector  responsivilitv  for  the  mission,  and  is  performed  as  closely  as 
practical  to  the  time  of,  and  under  the  same  conditions  as,  the  mission 
being  recorded.  This  calibration  is  verified  by  use  of  a  premission 
calibration. 

The  radiometric  calibration  is  performed  using  the  calibration 
collimator  designed  especially  for  this  purpose.  Basically,  the  colli¬ 
mator  consists  of  a  graybody  source  controllable  over  the  range  from 
about  300  to  1200  K,  a  set  of  apertures  that  may  be  placed  in  front  of 
Che  graybody,  and  a  system  of  collimating  optics  that  present  to  Che 
sensor  a  virtual  source  at  infinite  range.  While  the  GBM  Cracking 
pedestal  is  held  at  the  central  azimuth  location  and  the  GBM  primarv 
telescope  is  moved  in  elevation,  the  giaybody  image  is  slowly  stepped  up 
or  down  the  focal  plane  over  the  elevation  field  of  view. 

The  calibration  coefficients  have  been  derived  using  these  data 
and  various  assumptions.  Many  of  Che  assumptions  have  been  verified  to 
within  limits  by  analysis  of  more  extensive  calibrations,  such  as  tnose 
Chat  were  performed  during  laboratory  tests  prior  to  shipment  of  the  GBM 
system  to  Roi  Namur  for  installation. 

One  assumption  used  in  calibration  is  that  the  GBM  system  is  a 
linear  system,  since  evidence  accumulated  thus  far  suggests  chat  the 
sensor  response  is  nearly  linear  over  its  full  dynamic  range.  Another 
important  assumption  used  in  calibration  is  Chat  Che  response  of  the 
system  to  a  point  target  can  be  described  by  a  triplet  of  voltage 
values.  This  is  apparent  from  observations  that  Che  voltage  values  from 
a  point  target  are  frequently  detected  by  more  than  one  detector,  while 
hits  on  three  adjacent  detectors  are  very  comnon. 

B.2  THEORY 

Calibration  of  the  sensor  is  based  on  the  premise  that  one  can, 
through  a  combination  of  measurements,  calculate  the  effective 


B-1 


irradiance  present  at  the  sensor.  The  calibration  method  applied 
requires  only  that  the  energy  distributions  in  mission  and  calibration 
source  images  be  similar.  This  allows  data  obtained  with  both  focused 
and  moderately  defocused  systems  to  be  reduced.  The  method  allows  the 
total  irradiance  to  be  determined  even  when  a  portion  of  the  energy 
falls  on  an  inactive  detector  element. 

The  basic  theory  employed  in  computing  the  effective  irradiance 
produced  by  an  unresolved  source  is  based  on  the  condition  that  the 
source  irradiance  h^j  during  a  single  calibration  run  need  only  remain 
constant  during  the  recording  of  any  voltage  triplet.  These  voltage 
triplets  thus  satisfy  an  equation  of  the  form 

*bKo-i)i'  ’bji-  ''bO+Ui^  •  b.  ■  (B-I) 

where 

V.  . .  «  measured  waveform  voltage  above  the  local 
^  baseline  of  the  jth  detector  in  spectral 
band  b  of  the  1th  pulse 

h^  "  source  Irradiance  at  the  sensor  In  spectral 
band  b. 

This  fxinctlon^^  has,  assuming  only  linearity  of  detector  response,  the 
property 


J(f  C5V, 


b(j-l)’  °^b(j+l)^  ’  ^^b^\(j-] 


ba+D- 


where  5  is  any  arbitrary  factor.  Choosing  S  equal  to  the  square  root  of 
the  sum  of  the  squares  of  the  triplet  voltage  yields 


^'''b(J-l)'  ’bj’  ’'bU+l)^ 


where 


+  +  vf 


b(j+l)- 


“bj  -  ''b(j-i)/l''lbj 
\}  ■  ’bj'l’lbj  • 


’bU+D'I’lbJ 


This  choice  of  S  permits  to  be  expressed  as  the  product  of  two  factors 
one  of  which  depends  only  on  the  direction  cosines  in  the  voltage  domain. 
Expressed  in  this  manner  represents  a  universal  function  that  is  inde¬ 
pendent  of  the  particular  irradlance  used  in  obtaining  the  calibration 
data.  The  unknown  irradlance  can  therefore  be  obtained  by  multiplying 


b.3  ^ise  considerations 

Assuming  that  errors  in  the  measurements  caused  by  random  fluc¬ 
tuations  (noise)  are  additive  and  normally  distributed  about  the  mean 
value,  then  each  member  of  a  measured  voltage  triplet  may  be  represented 
in  the  form  V  =  V  +  Av,  where  V  is  the  true  (mean)  value  and  Av  is  the 
random  error,  which  is  normally  distributed  having  the  variance  £7^  . 
The  variance  may  be  different  for  each  detector  element. 

This  leads  to  the  concept  of  a  "resolution  cell"  in  the  voltage 
domain.  The  cell  is  ellipsoidal  in  shape  and  is  described  by  the  func¬ 
tion 


where  a*  ■  for  detector  J,  etc.  The  value  of  the  constant  in  Equation 
B-4  is  somewhat  arbitrary,  but  if  the  50*  point  is  selected,  the  value 
is  2.36606.  and  there  is  a  50Z  probability  that  the  observed  triplet 
V(j-l),  Vj .  V(j+i)  falls  within  the  ellipsoid  so  defined. 

The  existence  of  a  resolution  cell  of  finite  size  in  voltage 
space  suggests  that,  within  the  allowed  finite  ranges  of  voltages,  there 
are  only  a  finite  number  of  distinguishable  resolution  cells.  In  par¬ 
ticular,  if  the  number  of  voltage  triplets  measured  during  calibration 
exceeds  the  number  of  distinguishable  resolution  cells  in  voltage  space 
as  determined  by  random  noise,  it  is  justified  to  combine  some  of  the 
measured  triplets  by  averaging  to  reduce  the  total  number. 

Therefore,  for  each  pair  r,  s  (r  *  s,  s  ■  1  to  n,  where  n  is  the 
total  number  of  voltage  triplets  for  one  detector  triplet),  the  value  of 
Q«g  is  found  by 
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where  W  is  the  smaller  of  W_  and  W  .  The  inclusion  of  the  factor  W  in 
X  r  s  X 

the  formula  for  reflects  the  shrinking  size  of  the  resolution  cell  when 
points  are  combined  and  makes  it  less  likely  that  two  points  having  large 
and  nearly  equal  weights  will  be  combined.  The  general  effect  is  to  cause 
"clusters”  of  points  in  voltage  space  to  coalesce  to  a  single  point  each, 
with  the  weights  of  the  final  set  corresponding  to  the  number  of  points  in 
each  original  cluster.  Let  denote  the  smallest  calculated  above. 

If  ^  2.36606,  then  the  points  k  and  I  are  combined  by  averaging  accord¬ 
ing  to  the  formulas 


W  -  , 


^(j-1)  *  ^'^k  ^(j-l)k 


^J+l)k  ^ 


fB-b) 


The  two  original  points  k  and  I  are  then  replaced  by  the  single  average 
point,  thereby  reducing  the  total  number  by  one.  The  process  is  repeated 
until  >  2.36606. 


When  no  more  points  can  be  combined,  the  resultins  triplet 
voltages  are  converted  to  direction  cosines  -bjk’  *^bjk^  which 

together  with  the  voltages  constitute  the  calibration  table  for 

the  center  detector. 


B.4  CALIBRATION  FACTORS 

The  calibration  factors  were  derived  by  measuring  the  value 

of  the  waveform  voltage  above  the  local  baseline.  These  factors  were 
obtained  as  the  constant  image  of  the  collimating  graybody  source  slowlv 
moved  down  each  detector  in  elevation.  Different  responsivities  were 
measured  when  the  secondary  mirror  scanned  the  image  across  the  focal 
plane  right  to  left  and  when  it  scanned  left  to  right  (see  Figure  2-? 


for  focal  plane  detector  arrangement'.  This  direction-correlated  effect 
was  accounted  for  in  extraction  of  the  intensities. 

Normally,  mission  data  are  converted  to  engineering  units  on  the 
basis  of  calibration  data  recorded  immediately  before  and/or  after  the 
mission.  However,  in  the  case  of  mission  3TM-13W.  calibration  data  were 
recorded  only  on  two  detectors  of  the  LVIR  array  and  two  detectors  of 
the  SWIR  array.  The  reduction  process  requires  a  calibration  for  each 
detector  separately;  therefore,  the  calibration  data  recorded  for  mis¬ 
sion  GT  134M  were  substituted  in  the  conversion  of  measured  voltages  to 
units  of  irradiance.  This  calibration  file  was  selected  because  it  was 
of  good  quality  and  was  close  in  time  to  STM-1 3W. 

To  compensate  for  the  relative  responses  of  the  detectors 
between  the  calibrations,  it  was  necessary  to  develop  and  apply  a  set  of 
correction  factors.  These  procedures  are  outlined  in  Figure  B-1,  As 
shown  in  the  figure,  the  substituted  calibration  data  were  processed  in 
the  normal  manner  to  generate  a  calibration  file.  Voltage  peaks  for  the 
mission  calibrations  were  read  by  hand  from  strip  charts  and  converted 
to  irradiances  by  means  of  the  substituted  calibration  file.  Since  the 
effective  irradiances  during  the  mission  calibrations  are  known,  the 
departure  of  the  "calibration"  irradiances  from  this  value  represents 
the  error  to  be  corrected  using  suitable  correction  factors. 

Correction  factor®  C”  defined  by 

were  developed,  where  i  ■  1,  2,  ....  31  Identifies  the  particular  LWIR 
detector  and  ±  indicates  the  direction  of  scanning.  The  symbol 
denotes  the  average  value  of  the  "calibration"  irradiance  for  the  1th 
detector  and  the  ±  scan  direction.  Plnally,  the  correction  factors  are 
used  together  with  the  calibration  data  file  to  convert  the  mission 
voltages  to  irradiance  values. 
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AFPERDIX  C.  METEOROLOGY 


C.  1 


RADIOSONDE  BALLOON  FLIGHTS 


For  each  GBM  mission,  radiosonde  balloon  flights  were  conducted 
to  collect  meteorological  dat  supporting  the  Army  Optical  Station.  These 
data  included:  wind  direction,  wind  speed,  temperature,  dew  point, 
pressure,  density,  absolute  and  relative  humidity,  index  of  refraction, 
visibility,  and  wind  shear  as  a  function  of  altitude  between  12  and 
100,000  ft. 

Those  data  of  most  interest  in  the  analysis  of  clouds  are  the 
temperature  and  dew  point  profiles.  These  are  plotted  for  the  cases  of 
interest  between  0  and  40,000-ft  altitude  in  Figures  C-1  through  C-18. 


C.2 


CLOUD  RANGE  ESTIMATES 


A  rough  estimate  of  the  range  to  an  observed  cloud  can  be 
determined  using  the  boresight  tracking  data  of  Appendix  C  and  the 
Temperature  and  Dew  Point  Profiles.  These  estimates  are  determined 
assuming  that  the  temperature  and/or  dew  point  inversions  as  illustrated 
in  Figure  C-19  indicate  the  average  altitude  for  cloud  formations.  Then, 
using  the  boresight  elevation  of  the  cloud  measurement  and  simple 
trigonometry,  a  rough  estimate  of  the  cloud  range  is  made.  The  range 
error  associated  with  such  an  estimate  is  given  by, 


sin0 


TAN^e 


(C-1) 


where  the  elevation  angle  error  is  insignificant  compared  to  the 
altitude  error  ;  hence  Equation  C-1  reduces  to 


where 


sin0 


(C-2) 


altitude  error  (1  sigma) 
measured  elevation  angle  of  cloud. 
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This  appendix  contains  polar  plots  (Figures  D-I  through  D-8) 
o:  the  boresight  track,  and  sun  position  for  the  cases  of  interest.  The 
plot  center  represents  the  sensor  position  and  the  boresight  zenith. 
Each  bold  circle  represents  a  declination  of  10  deg  from  the  zenith, 
«itn  the  horizontal  represented  by  a  dashed  line.  The  sun  position  is 
also  snowr.  on  those  cases  for  which  it  appears  less  than  108  deg  from 
the  zenith. 
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